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Abstract 

Both the Dzyaloshiskii-Moriya interaction and the exchange-striction are shown to affect dynami- 
© : cally the magnetoelectric excitations in the perovskite mnltiferroic RMn0 3 . The exchange-striction 



results in a biquadratic interaction between the spins and the transverse phonons, giving rise to 
quantum fluctuations of the ferroelectric polarization P. This leads to low- lying phonon modes 
that are perpendicular to P and to the helical spins at small wave vector but are parallel to P at 
a wave vector close to the magnetic modulation vector. For spin-1/2 helimagnet, the local polar- 
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ization can be completely reversed by the spin fluctuation, and so does the direction of the on-site 
spin chirality, which allows for a finite differential scattering intensity of polarized neutrons from 

a 

a cycloidal magnet. 
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I. INTRODUCTION 

Details of the coupling mechanisms of the magnetic and the ferroelectric order in multifer- 
roics is currently under active research. This is due to the fundamental physics involved and 
to promising technological applications |l|. Our focus here is the perovskite multiferroics 
RMn03 with R = Tb, Dy, Gd and Eui-^Y^ that have incommensurate spiral spin structure 
2] . The experimental finding is that RM nO"3 has a helical magnetic order and a finite ferro- 
electric (FE) polarization. The driving mechanisms of this ordering is an interplay between 
the exchange interaction and the Dzyaloshiskii-Moriya (DM) interaction. Specifically, the 
spin-orbit coupling with a strength a related to the <i(p)-orbitals of the magnetic (oxygen) 



ions results in the FE polarization 



SEi 



4] P = acij x (Sj x Sj). Cij is a unit vector connecting 



the sites i and j. Generally, it is to be expected that the magnetoelectric (ME) coupling will 
affect not only the material static properties but also the dynamical response. Based on the 
spin-current model, the dynamical properties of DM interaction were studied in Ref.p-tZ]. 
A novel magnon-phonon excitations so-called electromagnon, was theoretically predicted. 
When the spiral plane rotates with respected to the axis of the helical wave vector, so does 



the the induced electric polarization, which couples the magnetic excitation to t 
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re electric 



field E of the radiation in the direction perpendicular to the spin spiral plane [5j|. Experi- 
mental observations in RMnOs 8| and Euo.75Yo.25Mn03 [9( seems to be consistent with this 
finding. However, a detailed study of the terahertz spectrum of Eui_ x Y x Mn03 [l0( revealed 
that infrared-absorption along the spontaneous polarization direction is also possible, which 
is not explained by theory mentioned above. This violation suggests that the static and the 



dynamic ME coupling may be different 11J. We carried out a detailed investigations of the 
dynamical properties of the multiferroics and find that both, the DM interaction and the 
(super) exchange striction play an essential role and need to be taken into account. 



II. THEORETICAL MODEL 



We consider a one-dimensional spin chain along the z-axis with a frustrated spin interac- 

nn 

tion. An effective model that captures the spin-phonon coupling [5|, H2[ has the Hamiltonian 



H = H s + H DM + H p (1) 

+ ^ J 2 (n - r m )3i ■ S m 

{I'Tfl/nnn 

H D m = -A J^ u, ■ [e z x (Si x S|+i)] 

i i 

where the notation (ij) nn indicates nearest-neighboring (nn) % and j, and (lm) nnn corre- 
sponds the next-nearest-neighboring (nnn) I and m. The competition between the nn 
ferromagnetic interaction (J x < 0) and the nnn antiferromagnetic interaction (J 2 > 0) 
leads to magnetic frustration and realizes a spiral spin ordering with the wave vector 



cosQ = — Ji/4J 2 [13Hl5|. H p describe optical phonons. The spin-phonon interaction Hdm 
originates from a spin-orbital coupling and breaks the inversion symmetry along the chain. 
Minimizing the energy yields the condition of the atomic displacement and the local spin- 
configuration, Uj = ^e z x (Si x Sj+i). Particularly if the zx helical spins along the chain, i.e. 
Sj = S(siniQ, 0, cosiQ), an uniform electric polarization P along the x direction is induced 
by the condensation of the transverse optical (TO) phonons, P = e*uo = — e*^f- sinQe^ 
with a Born charge e*. Generally, u x cannot be softened through the hybridization between 
the TO phonons and the magnons because of k/M ^> JS. The spontaneous FE polarization 
P x is frozen at — eu in the ferroelectric phase. However, after accounting for the superex- 
change striction, we have transverse acoustic (TA) phonons, which induces the fluctuation 
of the polarization hybridized with the spin bosons and soften thus the transverse phonon 
behavior. 

Considering small atomic displacements perpendicular to the chain, uj- ■ e z = 0, the 
exchange energy J falls off as a power law with the separation of the magnetic ions 

^{{v.-v^^J^l^-fini-nff] (2) 



where 7 is in the range of 6 — 14 16|. The emerging TA phonon mode is coupled to the spins 
with the bi-quadratic interaction — 71 ,2^1,2 ( u i~ ~ u f) 2 {Si ' S,-). This dynamical coupling does 
not contribute any additional static electric polarization but induces the fluctuation of the 



electric dipole moment due to the low frequency excitation modes of TA phonon. We write 
explicitly the atomic displacements into two parts: (i) the statical part u« = (uq, 0, 0) and 
(ii) the dynamical part #Uj = (—Suf, Suf, 0). Retaining terms up to the second order in the 
quantum fluctuation, the spin-current model delivers the following coupling terms 

H DM = -XS cos Q J2 K(Si+i - Sf) 

i 

-\S 22 $ uV i(Si cos Qi+l ~ &i+l cos Qi) ( 3 ) 

/ 

in the rotated spin frame: Sf = SfcosiQ + SfsiniQ, Sf = Sf, and Sf = —SfsiniQ + 
Sf cos iQ. 

ILL RESULTS AND ANALYSIS 

In spin-1/2 multiferroics, such as LiCu202 ll4j, the spin fluctuations may spontaneously 
reverse the local spin. Defining the vector of spin chirality as the average of the outer 
product of two adjacent spins c« = (s$ x Sj +1 )/|s.; x s.; +1 |, in the RMn0 3 -type multiferroics the 
direction of local FE polarization is determined by the on-site spin chirality. The dynamical 
DM interaction Eg. ([5]) yields the coupling term between the spin and the spin- chirality in 
the spin-1/2 multiferroics, £\ c?(lf +1 — sf) = J2 i sf{cf_ t — cf ), which indicates that when the 
spin at site i is flipped, Si — > — §i, the direction of spin- chirality q and q_i are also reversed. 
Assuming all spins point along their corresponding classical directions in the ground state 
of the spin-1/2 helical magnet as in NaCu202, where a J\ — J2 spin model provides a good 
description of the helix state [15| . So the spin interaction can be ferromagnetically given as 
—J s (Q)si ■ Sj where Q is taken as the pitch angle along the chain. An effective model that 
describes the interplay between the helical spin and spin-chirality has the form 

H sc = - X)Wi • % + J c c t ■ dj) - 7 Yl § i( d i-i ~ %)■ ( 4 ) 

i,j i 

The Hilbert space can be considered as the tensor product space \i) — > \sf) s <8> \cf) c . Now 
if the spin at site % is flipped, the spin and spin-chirality excitations are mixed due to the 
spin-phonon coupling. The expected value of spin-chirality is given by 

(c) = l-(s), (5) 



which is less than one. The experimental data for a finite differential scattering intensity 
of polarized neutrons from LiCu202 [14J suggests (c) ~ 0.3 which is consistent with the 
estimated value (c) = 0.44 based on the ordered moment, 0.56/i^ per magnetic copper site 

m 

For magnet RMnOs, the helical spin ordering occurs, corresponding to the condensation 
of the spin bosons. By using the standard linear-spin-wave approximation, a dynamical 
magon-phonon interaction reads, 

H D m = -AS 1 cos Q ^2 ?>UqS%(cos q-1) 

<? 

-xsj2^ y A Q (^ iQ -^ q±Q) )/^ (6) 

q 

8u y is hybridized with the spin at q ± Q (optical magnons), but Su x is coupled to S x at q 
(acoustical magnons). The polarization correlation functions are given as 



< 5u x q \5u x > 



u 2 -cj 2 



M[w*-w>»+w?)+a;>J-fc;y]' 
« 8u\\8u\ »= — - A2g 3 —-■ 

where u p is the frequency for the transverse phonon, u s (q) is the energy dispersion of the 
spin-excitation, u' sp (q) = [2(A(q) - 2B{q)) (\ 2 S 3 cos 2 Q(l - cosg))/^] 1 / 2 , and G s (q ± Q) = 
(A(q ±Q) + 2B(q ± Q))(l - cos(g ± 2Q)/(co 2 - co s (q ± Q)) with 

A(q) = — JifcosQ + -(1 + cosQ) cosg] 

- J 2 [cos2Q + -(l + cos2Q)cos2g] 

A^sin^Q 
+ ^fe (2-cosg) (7) 

B(q) = -r(l — cosQ) cosq H (1 — cos2Q) cos2g 

A 2 S 2 sin 2 Q 
^ cos ^ ( 8 ) 

At small wave vectors, q ~ and u p ~ \/k/M, the TA phonon is decoupled from spins. The 
antisymmetric DM interaction dominates over the spin-phonon coupling. 5u y Q is coupled via 
(Sq — Sq) to the rotation of the spin plane and the direction of the polarization along the 
chain. However, at a wave vector close to the magnetic modulation vector, i.e. q ~ Q, both 
the symmetric and antisymmetric magnetoelectric interaction respond to the fluctuations 



of the polarization. Especially, in the direction parallel to the FE polarization P, there is a 
low frequency range around bj x i = u s (Q) where u x couples resonantly to light. Introducing 
an easy-plane spin anisotropy D(S y ) 2 into the spin system, we observe nearly the same 
low-frequency behavior of the polarization correlation functions oj x i ~ V ' JSD m u y _. These 
conclusions are also qualitatively consistent with experiment observations for Eui_ x Y x Mn03 



io|. 



IV. SUMMARY 

In conclusion, we studied the origin of the magnetoelectric dynamics in the orthorhombic 
perovskite RMnOa. At a small wave vector, the DM interaction determines the low-frequency 
behavior of the phonons. For a wave vector close to that of the magnetically modulated 
structure, the exchange striction induces fluctuations in the FE polarization, and additional 
low-lying mode parallel to the FE polarization emerges. Due to the dynamical Dzyaloshiskii- 
Moriya interaction, the spin-chirality is strongly coupled to the spin fluctuation which implies 
a large quantum ffucatuation of the spin-chirality in the ordered spin- 1/2 system and results 
in a finite scattering intensity of polarized neutrons from a cycloidal helimagnet. 
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